Microinstabilities at perpendicular collisionless shocks: A comparison of full particle 
simulations with different ion to electron mass ratio 
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A full particle simulation study is carried out for studying microinstabilities generated at the 
shock front of perpendicular collisionless shocks. The structure and dynamics of shock waves are 
determined by Alfven Mach number and plasma beta, while microinstabilities are controlled by the 
ratio of the upstream bulk velocity to the electron thermal velocity and the plasma-to-cyclotron 
frequency. Thus, growth rates of microinstabilities are changed by the ion-to-electron mass ratio, 
even with the same Mach number and plasma beta. The present two-dimensional simulations show 
that the electron cyclotron drift instability is dominant for a lower mass ratio, and electrostatic 
electron cyclotron harmonic waves are excited. For a higher mass ratio, the modified two-stream 
instability is dominant and oblique electromagnetic whistler waves are excited, which can affect the 
structure and dynamics of collisionless shocks by modifying shock magnetic fields. 
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PACS numbers: 52.35.Tc; 52.35.Qz; 52.65.-y; 52.65.Rr; 



I. INTRODUCTION 

Collisionless shocks have been investigated by full par- 
ticle (Particle-In-Cell) simulations for more than four 
decades since early 1970's The full particle method 
handles both electron-scale microphysics and ion-scale 
shock non-stationarity simultaneously because both elec- 
trons and ions are treated as individual charged particles. 
However, a reduced ion-to-electron mass ratio {mi/me) 
is commonly used in the full particle method to reduce 
the computational cost. 

It is well known that structures and dynamics of shock 
waves are determined by the following two dimensionless 
parameters, i.e., Alfven Mach number 
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and the ratio of the thermal plasma pressure to the mag- 
netic pressure (plasma beta) 



/3e = 
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where c, u, ujp, lOc, vt and m represent the speed of light, 
upstream bulk velocity, plasma frequency, cyclotron fre- 
quency, thermal velocity and mass, respectively, with the 
subscripts "i" and "e" being ion and electron, respec- 
tively. Here the subscript "1" denotes upstream. 
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It is well known that a part of incoming ions are re- 
flected at the shock front of perpendicular and quasi- 
perpendicular shocks. The reflection of ions results in 
the deceleration of incoming electrons so that the con- 
servation of the total current (the zero current condition 
in the shock normal direction) is satisfied. Consequently, 
a relative bulk velocity between the incoming electrons 
and the incoming/reflected ions arises. Previous self- 
consistent kinetic simulation studies revealed that there 
exist various types of microinstabilities at the shock front 
of perpendicular and quasi-perpendicular shocks in which 
reflected ions play essential roles in dissipative processes. 
The free energy source of these microinstabilities is the 
relative drift (current) between incoming electrons and 
incoming/reflected ions. 

The ratio of the relative drift velocity to the electron 
thermal velocity is important for controlling the type of 
microinstabilities (or the effect of electron thermal damp- 
ing to waves). This ratio is proportional to the ratio of 
the upstream bulk velocity to the electron thermal ve- 
locity. From the equations of Alfven Mach number and 
electron beta, we obtain 



(3) 



This relation means that the ratio of the upstream bulk 
velocity to the electron thermal velocity becomes larger 
with larger Alfven Mach number, smaller electron beta, 
or smaller mass ratio. Note that the actual amount 
of free energy relative to the electron thermal energy 
(miU^i/ nieVf^i) is independent of the mass ratio. 

If the relative velocity between electrons and ions ex- 
ceeds the electron thermal velocity, electrostatic waves 
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are excited by current-driven instabilities such as the 
Buneman-type instability (BI) [21 or the electron cy- 
clotron drift instability (ECDI) [3, i]. At high-Mach- 
number perpendicular shocks, the relative velocity be- 
tween incoming electrons and reflected ions commonly 
becomes much faster than the electron thermal velocity. 
Then, the BI becomes dominant, and electrostatic waves 
are excited at the upper hybrid resonance frequency 
At lower- Mach- number {Ma < 10) perpendicular 
shocks, the relative velocity between incoming electrons 
and incoming/reflected ions becomes close to the elec- 
tron thermal velocity. Then, the growth rate of the BI 
becomes small because of damping by thermal electrons, 
and the ECDI becomes dominant, which excites electro- 
static waves at multiple electron cyclotron harmonic fre- 
quencies Q. When the relative velocity between incom- 
ing electrons and incoming/reflected ions becomes slower 
than the electron thermal velocity at lower-Mach-number 
perpendicular shocks, high-frequency electrostatic waves 
are not excited due to damping by thermal electrons, 
and the modified two-stream instability (MTSI) becomes 
dominant 0,i,i,liSlil|. Then, obliquely propagating 
electromagnetic whistler mode waves are excited at a fre- 
quency between the electron cyclotron frequency and the 
lower hybrid resonance frequency [H, [H, [ij, [H, [13, 

m. 

The previous theoretical work has shown that the 
growth rate of the MTSI is strongly affected by the ion- 
to-electron mass ratio (mi/TOg), but not by the ratio 
of the electron plasma to cyclotron frequency (wpe/wce) 
(l^ . It has also been confirmed by the previous one- 
dimensional (ID) full particle simulation studies that the 
mass ratio affect the shock dynamics as well as the mi- 
croinstability jisl [Tsj . In the present study, we vary 
the mass ratio and frequency ratio as shown in Table 
I. Even when the Alfven Mach number and the electron 
beta are fixed, the ratio of the plasma frequency to cy- 
clotron frequency (wp/wc) can take an arbitrary value, 
because c/vte or ui/c can also take arbitrary values. The 
frequency ratio may also control the type of microinsta- 
bilities by determining the linear dispersion relation of 
background plasma. It should be noted that the effect 
of the frequency ratio is not included in hybrid simula- 
tions where ions are treated as particles while electrons 
are treated as a (massless) fluid. That is, the ion inertial 
length {c/ujpi) and the ion gyro radii (vu/ujci and ui/cud) 
are defined but ujpi/uJci oo and c/vu — oo are assumed 
in hybrid simulations. Thus the full particle simulation 
is the unique approach to study the effect of both mass 
ratio and frequency ratio. 

In the present study, we perform a series of two- 
dimensional (2D) full particle simulations of a perpen- 
dicular (the shock normal angle 9b„ — 90°) coUisionless 
shock with a low Mach number {Ma = 6) and a mod- 
erate beta {/3i = /3e — 0.32). We examine the impact of 
the mass ratio on microinstabilities at the shock front by 
fixing either Upn/tUcii or ojpei/wcei as shown in Table I. 



II. SIMULATION SETUP 

We use a 2D electromagnetic full particle code in which 
the full set of Maxwell's equations and the relativistic 
equation of motion for individual electrons and ions are 
solved in a self-consistent manner. The continuity equa- 
tion for charge is also solved to compute the exact cur- 
rent density given by the motion of charged particles 
[l9| . In the present simulation, the simulation domain 
is taken in the shock-rest frame [l^, [2l|, [g^l . The shock- 
rest-frame model is achieved by using the "relaxation 
method" @, H^, H^l , in which a coUisionless shock is ex- 
cited by an interaction between a supersonic plasma flow 
and a subsonic plasma flow moving in the same direction. 

The initial state consists of two uniform regions sep- 
arated by a discontinuity. In the upstream region that 
is taken in the left-hand side of the simulation domain, 
electrons and ions are distributed uniformly in space 
and are given random velocities {vx,Vy,Vz) to approxi- 
mate shifted Maxwellian momentum distributions with 
the drift velocity ui, number density rii = eoWeWpgi/e^, 
isotropic temperatures T^i = niev'^gi and Tn = rriiV^^^, 
where e is the charge. Here, the subscripts "1" and 
"2" denote "upstream" and "downstream", respectively. 
The upstream magnetic field Bqi with a magnitude of 
TOeWcei/e is also assumed to be uniform. The down- 
stream region taken in the right-hand side of the simula- 
tion domain is prepared similarly with the drift velocity 
U2, density n2, isotropic temperatures Te2 and Ti2, and 
magnetic field i?02- 

We take the simulation domain in the x-y plane and as- 
sume an in-plane shock magnetic field {Byo). The shock- 
normal magnetic field B^o is set to be zero. That is, 
the shock-normal angle 6'b„, which is defined as the an- 
gle between upstream magnetic field and shock-normal 
direction, is set as 90° in the present study. As a 
motional electric field, a uniform external electric field 
Ezo = UxiByoi = Ux2By02 is applied in both upstream 
and downstream regions, so that both electrons and ions 
drift along the x axis. At the left boundary of the simu- 
lation domain in the x direction, we inject plasmas with 
the same quantities as those in the upstream region, while 
plasmas with the same quantities as those in the down- 
stream region are also injected from the right boundary 
in the x direction. We use absorbing boundaries to sup- 
press non-physical reflection of electromagnetic waves at 
both ends of the simulation domain in the x direction 
[2^, while the periodic boundaries are imposed in the y 
direction. 

In the relaxation method, the initial condition is given 
by solving the shock jump conditions (Rankine-Hugoniot 
conditions) for a magnetized two-fluid isotropic plasma 
consisting of electrons and ions [26]. In order to de- 
termine a unique initial downstream state, we need 
given upstream quantities u^i, tOpei, "^cei, Vtei, vui and 

Ti2/Te2- 

In the present study, we fix the Mach number and 
plasma beta to Ma = 6 and fie = Pi = 0.32. As shown in 
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TABLE I: Parameters for difTerent simulation runs. 



Run 


mi/ me 


CJpel/I^cel 


I^pil/I^cil 


Uxl/Vtel 


C/Vtel 


7 / \ 

til/ AOel 


A" 


625 


4.0 


100.0 


0.6 


10.0 


250 


B 


256 


4.0 


64.0 


0.9375 


10.0 


160 


C 


100 


4.0 


40.0 


1.5 


10.0 


100 


D 


25 


4.0 


20.0 


3.0 


10.0 


50 


E 


256 


6.25 


100.0 


0.9375 


15.625 


250 


F 


100 


10.0 


100.0 


1.5 


25.0 


250 


G 


25 


20.0 


100.0 


3.0 


50.0 


250 
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Table I, we perform 7 different runs by varying the mass 
ratio {rrii/me) and frequency ratio. We fix tOpe/uJce in 
Runs B-D. Then, the electron parameters are fixed while 
the ion parameters are changed. On the other hand, we 
fix idpi/cjci in Runs E-G. Then, the ion parameters are 
fixed while the electron parameters are changed. 

The size of the simulation domain is Lx — 32Zii and 
Ly = hi iu the shock-normal and shock-tangential direc- 
tions, respectively, where In = c/tUpn is the ion inertial 
length. The grid spacing is set to be Ax = Ay = Xoei 
where XdcI is the electron Debye length upstream. On 
the other hand, the time step is set to be cAt/Ax = 0.5. 
These mean that the computational resource (the size 
of the simulation domain and the number of time step) 
increases with the mass ratio when the electron param- 
eters are fixed (Runs A-D), while the computational re- 
source becomes constant when the ion parameters are 
fixed (Runs E-G and A). We used 16 pairs of electrons 
and ions per cell in the upstream region and 64 pairs 
of electrons and ions per cell in the downstream region, 
respectively, at the initial state. 



III. RESULTS 

A. Simulation Results 

Figure 1 shows the evolution of the perpendicular 
shock in Runs A-G. The snapshots of the magnetic field 
Bx, By, Bz components, the electric field component, 
and the corresponding x — phase-space distribution 
function of ions at uJdit = 10.2 in Runs A-G are plotted 
in Figures A-G, respectively. The position is normalized 
by the ion inertial length In = c/ujpn. The magnitudes 
of magnetic field and electric field are normalized by the 
initial upstream magnetic field Byoi ^^nd the motional 
electric field E^q, respectively. In the panel (a'), the tan- 
gential component of the magnetic field By in Run A 
is plotted as a function of position x and time t. Note 
that the tangential component of the magnetic field By in 
Runs B-G shows almost the same evolution as the panel 
(a'). 

In the present shock-rest-frame model, a shock wave 
is excited by the relaxation of the two plasmas with dif- 
ferent quantities. Since the initial state is given by the 
shock jump conditions for a "two-fiuid" plasma consist- 



ing of electrons and ions, the kinetic effect is excluded in 
the initial state. Thus the excited shock slightly moves 
upstream. The shock front appears and disappears on a 
timescale of the ion cyclotron period in Figure la', which 
corresponds to the cyclic shock reformation P, [23, • 

In all the runs, a certain part of the incoming ions is 
reflected upstream at the shock ramp and the shock foot 
region is formed, which is called the broadening phase 
of the reformation. Consequently, microinstabilities are 
generated due to the relative drift between the reflected 
ions and the incoming electrons, and ion phase-space vor- 
tices in association with wave generation are formed as 
seen in the ion x — phase-space plots. The spatial 
size of the foot region corresponds to the gyro radius of 
reflected ions, which seems to be longer with respect to 
the ion inertial length as the mass ratio becomes smaller 
when the electron parameters are fixed (Runs A-D). It is 
found that the wavelength of the magnetic field B^ com- 
ponent in the shock normal direction becomes longer as 
the mass ratio becomes smaller. The excited waves in E^ 
and Bz components have the same wavelength in Runs 
A, B, C, E, and F, indicating the excitation of electro- 
magnetic waves. Also, the wavelength is shorter with the 
fixed LOpi/uJci than with the fixed Wpe/'^ce- However, the 
saturation amplitude of these waves is almost indepen- 
dent of the mass ratio (typically SB^ ^ 0.3i?o). On the 
other hand, the excited waves in Runs D and G have a 
shorter wavelength in E^ component than in B^ compo- 
nent, indicating the excitation of electrostatic waves in 
the electric field E^ component. The typical timescale 
of the saturation of instabilities is ~ 0.6/ujcii from the 
beginning of the broadening phase (ion reflection). Note 
that the typical timescale of the reformation is ^ 1.6/wcii 
in the present case. 

In Figures 2 and 3, we show spectra of Bz and E^ 
components, respectively, for Runs A-G. These spectra 
are obtained by Fourier transformation in the range of 
-6.5 < x/hi < -5.5, < y/l,i < 1, and 9.8 < u}„it < 
10.6. The intensity of Bz and E^ is normalized by Byoi 
and Ezo, respectively. In the top panels, the uj — — ky 
spectrum is reduced to the kx — ky space by integrating 
over uj. In the bottom panels, the uj — kx — ky reduced to 
the UJ — kx space integrating over ky. 

The spectra of both Bz and Ex show that there are 
several wave modes. An ion-scale component exists at 
a; ~ 0, fca; ~ and ky ~ 0, which may correspond to 
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L-mode ion-cyclotron waves due to the ion temperature 
anisotropy. However, this wave mode is not focused on 
in the present study because the present Fourier analysis 
cannot resolve ion scales. 

In all the simulation runs, waves with an electron-scale 
wavelength are excited almost in the shock-normal direc- 
tion {ky ~ 0). The excited waves propagate mostly in 
both (ice) directions, but the waves propagate in the —x 
direction are dominant. The phase velocity of the excited 
waves in the —x direction estimated from the uj — k^ spec- 
tra corresponds to the drift velocity of the reflected ions. 
The phase velocity of the excited waves in the +x di- 
rection corresponds to the drift velocity of the incoming 
ions. 

In all the runs, the frequency of the Bz component 
is lower than the upstream electron cyclotron frequency 
LOcei but is higher than the upstream lower hybrid reso- 
nance frequency {lulhri/^ccI ^ 0.04 in Run A, 0.06 in 
Runs B and E, 0.1 in Runs C and F, and 0.2 in Runs 
D and G). These properties agree well with the scenario 
that whistler mode waves propagating in the direction 
oblique to the upstream magnetic field are excited by 
the MTSI due to ion beams (incoming or reflected ions) 
across the shock magnetic field. In Runs D, F, and G, 
there are also enhancements of the component at fre- 
quencies higher than the upstream electron cyclotron fre- 
quency ujcei but is lower than the upstream upper hybrid 
resonance frequency {uJuHRi/'^pei = 4.1 in Run D, 10.0 
in Run F, and 20.0 in Run G). These properties agree 
with the scenario that electron cyclotron harmonic waves 
propagating in the direction almost perpendicular to the 
upstream magnetic field excited by the ECDI due to ion 
beams (incoming or reflected ions) across the shock mag- 
netic field. These results suggest that the ECDI is gen- 
erated when the ion-to-electron mass ratio is relatively 
small. 



B. Comparison with linear analysis 

A previous linear analysis and numerical simulation by 
Matsukiyo and Scholer [1^ suggested that there are two- 
types of instabilities due to both incoming and refiected 
ions. Hereafter, we refer to an instability due to incoming 
ions as "instability-1" , and refer to an instability due to 
refiected ions as "instability- 2" . 

In order to analyze the linear dispersion relation at 
the shock foot, we approximate the velocity distribution 
function with Maxwellian velocity distributions, that is 
given by 



exp 



2V? 



(4) 



For simplicity, we focus on the velocity distribution func- 
tion only in the shock- normal direction (x component). 
The velocity distribution function f(vx) is modeled by 
the three-component plasma consisting of incoming ions 
(il), reflected ions (12) and bulk electrons (e). Figure 4 



shows the typical velocity distribution functions of ions 
and electrons at the shock foot just before the genera- 
tion of instabilities in Runs A-G. Table 2 shows detailed 
parameters for the Maxwellian velocity distribution func- 
tions estimated from Fig. 4. The velocity, density and fre- 
quency are normalized by the upstream electron thermal 
velocity Vtei , the upstream density noi , and the upstream 
electron plasma frequency ujpei, respectively. 

The actual dispersion relation is solved by assuming 
that the ions are unmagnetized and the electrons are 
magnetized as in the previous studies [l3, [3 • 1^ should 
be noted that although both ion and electron components 
have a certain shift in the Vz space (not shown) , the effect 
of the drift velocity in the z direction is not included in 
the present analysis. 

Figure 5 denotes frequencies (upper panel) and growth 
rates (lower panel) of the instability-1 (fc > 0) and 
instability-2 (fc < 0) below the electron cyclotron fre- 
quency for a variety of wave propagation angles with 
respect to the ambient magnetic field, Oskj 3.S a func- 
tion of kx — ksmOs^. and ky = kcosOB^- For the 
above parameters, both of the two instabilities are un- 
stable at the phase velocities corresponding to the drift 
velocities of incoming and refiected ions. In all the runs, 
obliquely propagating waves with a frequency a; < ujcei 
are unstable, which corresponds to the MTSTl due to 
incoming ions and the MTST2 due to refiected ions. We 
found that the MTSIs have a maximum growth rate at 
kyki ^ 4. Although the linear growth rates of both 
MTSI-1 and MTSI-2 are of the same order in the linear 
analysis, the MTSI-2 is dominant in all the simulation 
runs. One major reason for this is that the wavenum- 
ber resolution in the shock tangential (y) direction is 
Aky = 6.28/;,i - 0.25/^61 in Run A, - 0.39/Zei in Runs 
B and E, - 0.63/Zei in Runs C and F, and - 1.26/Zei 
in Runs D and G, so that the simulation system allows 
the excitation of waves only with a larger wavenum- 



ber {ky 



kx/ tanOBf. > Aky) [18l |. The maximum 



growth rate occurs at uj/ujcei 0.14, kxlei ~ 1-26 in 
Run A, uj/ujcei ^ 0.26, k^lei ^ 1-57 in Runs B and E, 
uj/ujcei ^ 0.47, kxlei ~ 1-98 in Runs C and F, which 
are almost in agreement with the simulation result. For 
Runs D and G, the system does not have enough fre- 
quency and wavenumber resolutions to directly compare 
the theory and simulation. The result suggests that the 
MTSI-2 is generated for a wide range of the mass ratio, 
although the growth rate of the MTSIs becomes smaller 
as the mass ratio becomes smaller. 

Figure 6 denotes frequencies (upper panel) and growth 
rates (lower panel) of the instability-1 {kx > 0) and the 
instability-2 {k^ < 0) above the electron cyclotron fre- 
quency as a function of k^ for Runs C, D, F, and G. For 
the above parameters, both of the two instabilities are 
unstable at the phase velocities corresponding to the drift 
velocities of incoming and refiected ions. In theses runs, 
waves propagating in the direction perpendicular to the 
shock magnetic field are unstable at w '-^ nujcei, which 
corresponds to the ECDI-1 due to incoming ions and 
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TABLE II: Parameters for the Maxwellian velocity distribution functions estimated from Fig. 4. The velocity, density and 
frequency are normalized by the upstream electron thermal velocity Vtei, the upstream density tiqi, and the upstream electron 
plasma frequency ujpei , respectively. 



Run 




Vtii 


nil 


Vdi2 


Vu2 


ni2 




Vte 


He 






A 


0.42 


0.06 


1.68 


-0.32 


0.05 


1.38 


0.08 


1.94 


3.06 


1.75 


0.92 


B 


0.60 


0.11 


1.68 


-0.52 


0.07 


1.04 


0.17 


1.76 


2.72 


1.65 


0.92 


C 


0.96 


0.17 


1.71 


-0.87 


0.11 


1.27 


0.22 


1.80 


2.98 


1.73 


0.92 


D 


2.12 


0.28 


1.70 


-1.85 


0.21 


1.41 


0.32 


1.75 


3.11 


1.76 


0.92 


E 


0.58 


0.10 


1.72 


-0.52 


0.06 


1.10 


0.15 


1.89 


2.82 


1.68 


0.59 


F 


0.99 


0.15 


1.68 


-0.81 


0.10 


1.37 


0.21 


1.91 


3.05 


1.75 


0.37 


G 


2.15 


0.30 


1.68 


-1.42 


0.30 


1.40 


0.56 


1.86 


3.08 


1.75 


0.18 



the ECDI-2 due to reflected ions. The maximum growth 
rate of the ECDI-2 occurs at w ^ 3.4wcei, kxlei ^ —6.3 
in Run D, cj ^ B.lcjcei, k^lei ^ —12.5 in Run F, and 
Lo 3.2ujcei, k^lei —8.7 in Run G, which are almost 
in agreement with the simulation result. In Run C, on 
the other hand, the maximum growth rate of the ECDI-2 
occurs at ~ S.lojcei, k^lei ~ —10.6 in the linear anal- 
ysis, while there is no wave excitation by ECDI-2 in the 
simulation. This is because Run C does not have enough 
resolution for waves at = 10 /lei ~ l/Auei- We per- 
formed a high-resolution {Ax = O.lA^jei) ID simulation, 
and confirmed that the ECDI-2 is weakly generated. It 
should be noted that the linear analysis for Runs A, B 
and E shows a positive but small {-^/iOcei < 5) growth 
rate for the ECDIs at a wavelength shorter than the elec- 
tron Debye length {kxXoei > !)■ However, full particle 
simulations do not show the generation of the ECDIs 
even with a higher resolution. 

Although both ECDI-1 and ECDI-2 have a positive 
growth rate in the linear analysis, the ECDI-1 is not gen- 
erated in these simulation runs. Note that the simulation 
system has enough spatial resolution for ECDI-1 in Runs 
D, F and G (see Table 1). Here we present several reasons 
for this by using the temporal development of the foot 
region shown in Fig. 7. First, we analyze the velocity dis- 
tribution functions near the ramp-side of the foot region, 
while the ECDI-1 is generated from the upstream-side 
of the foot region. Second, waves excited in the ramp- 
side of the foot region quickly reach the shock front and 
the ECDI-1 does not have enough growth time. Third, 
the relative drift between the incoming ions and electrons 
arises later in the broadening phase. The incoming ions 
move through a distance of 2.4/ii in a time of OA/uJch 
(because of Ma = 6). As seen in Fig. 7, the structure 
of the foot region quickly changes on the timescale of 
0.1/iOcii, and the relative drift between the incoming ions 
and electrons exists in a narrow region of < 2lii. 



IV. CONCLUSION AND DISCUSSION 

The present 2D full particle simulations suggested that 
the ion-to-electron mass ratio affects microinstabilities in 
the foot region of perpendicular coUisionless shocks. Both 
simulation result and linear analysis clearly show that 



electromagnetic oblique whistler waves due to the MTSI 
are dominant with a large mass ratio, while electrostatic 
ECH waves due to the ECDI are dominant with a small 
mass ratio. 

The electrostatic ECH waves contribute to the genera- 
tion of nonthermal population of electrons [221 . The elec- 
tromagnetic oblique whistler waves contribute to the elec- 
tron heating in the direction parallel to magnetic fields, 
which slightly modifies the shock jump condition. The 
influence of the excited waves to ion heating is small in 
the present study. 

It should be noted that the MTSIs are generated even 
with small mass ratios (Runs B-G), which is seemingly in- 
consistent with the previous expectation |1J]. However, 
the previous linear analysis by Matsukiyo and Scholer 
[TJ] was done by using local (foot) plasma parameters as 
normalization factors, although the linear analysis here 
utilizes the upstream plasma parameters as normaliza- 
tion factors. The period of reformation is of the order 
of upstream ion gyro period which is often about three 
times longer than a local ion gyro period estimated by 
using the local magnetic field value in the foot. 

In most of the cases {mi /me > 100: except for Runs 
D and G), the MTSIs are predominant in terms of wave 
power compared with the ECDIs. As pointed out in the 
previous studies [3l|, HI] , ECDI saturates at a low level. 
Furthermore, it is confirmed in a two dimensional pe- 
riodic (local) simulation that the MTSI finally survives 
even though the ECDI has a larger growth rate [l^ if 
MTSI and ECDI are simultaneously present. Here, we 
confirmed such properties in our self-consistent 2D shock 
simulations. 

In the shock-tangential direction, large-amplitude fluc- 
tuations on the ion inertial scale are commonly exist at 
the shock front, which are known as the "ripples" (33| . 
Recent full particle simulations with a large system size 
suggested that shock-front ripples change the structure of 
shocks and enhance microinstabilities at the shock front 
[isl [22I . [29I Isoj . In the present study, however, the gen- 
eration of the ripples is not included because the system 
length in the shock-tangential direction is too short. The 
present linear analysis shows that the MTSIs have the 
maximum growth rate at k^^ln 4, which suggests that 
the oblique whistler waves can couple with L-mode ion- 
cyclotron waves at the shock overshoot that are a seed 



6 



perturbation of the ripples. A large-scale 2D simulation 
including shock-front ripples is ciirrently undertaken to 
understand the competition between the ripples and mi- 
croinstabilities, and is reported in a future paper. 
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FIG. 1: The evolution of the perpendicular shock in Runs A-G. (a-g) Magnetic field B^, By, B^ components, electric field 
Ex component, and corresponding x ~ Vx phase-space distribution function of ions at ujciit = 10.2 in Runs A-G. The solid line 
in the phase-space plot indicates the electron bulk velocity, (a') Tangential magnetic field By averaged over y as a function 
of position x and time t in Run A. The position and time are normalized by Ui and I/luch, respectively. The magnitude of 
magnetic field is normalized by the initial upstream magnetic field Byoi. The magnitude of electric field is normalized by the 
motional electric field Ezo. 
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FIG. 2: Spectra of Bz component for Runs A-G. The magnitude is plotted in a logarithmic scale from 10 ^Bj/oi to 10 ^Bj,oi- 



9 



(a) 'H,/fiic=G25 



(b) /?I.:/V,-256 

<jJp,i ■■(.■),, =4 
",ii/^Vii-&4 



(C) fl;,: -H~100 (d) /f7j/>T^^-25 



-1 -s -1 fi 



4 J 5 



n ii 

-; -1 -S -1 fi ■ 2 i A 5 



-it -■: - J 



2C- 



.1£ '■H I 



-i: - t -IS ■: & j 



ii 



1 

— J ? 

u 

5 



■ :■ ^ 




— -I 

1' 



I 



I II 



■■■ I II 



d 



(e) fl^f/''^-25e (f) ynj/>7?,,-100 

■■'f'>..j-10D Oi, i/t\ii-100 



(g) - 25 

Cl>p._/D},.|=100 



9.S<C0,i,/^ 10.6 ^ 
-6 5^ X//.. <-5.5 




FIG. 3: Spectra of Ex component for Runs A-G with the same format as Fig.2. The magnitude is plotted in a logarithmic 
scale from lO'^Ezo to lO'^^zO- 
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FIG. 4: Reduced velocity distribution functions of ions and electrons {fi{vx) and fe{vx)) averaged over y at the shock foot 
just before the generation of instabilities in Runs A-G. 
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FIG. 5: Linear dispersion relations of the MTSI-1 (fe > 0) and MTSI-2 [kx < 0) for Runs A-G based on the velocity 
distribution functions in Fig.4. Frequencies and growth rates are plotted for a variety of wave propagation angles with respect 
to the ambient magnetic field, 9bi,, as a function of kx = fcsin^sj. and ky = k cos 6b^- 
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FIG. 6: Linear dispersion relations of the ECDI-1 (kx > 0) 
and ECDl-2 {kx < 0) for Runs C, D, F and G based on the 
velocity distribution functions in Fig.4. 




FIG. 7: The x ~ Vx phase-space distribution function of ions 
at different times in Run G. The solid line in the phase-space 
plot indicates the electron bulk velocity. 



